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Abstract. We present self-consistent dynamical models for dust-driven winds of carbon-rich AGB stars. The 
models are based on the coupled system of frequency- dependent radiation hydrodynamics and time-dependent 
dust formation. We investigate in detail how the wind properties of the models are influenced by the micro-physical 
properties of the dust grains that are required by the description of grain formation. The choice of dust parameters 
is significant for the derived outflow velocities, the degrees of condensation and the resulting mass-loss rates of 
the models. In the transition region between models with and without mass-loss the choice of micro-physical 
parameters turns out to be very significant for whether a particular set of stellar parameters will give rise to a 
dust-driven mass- loss or not. We also calculate near-infrared colors to test how the dust parameters influence the 
observable properties of the models, however, at this point we do not attempt to fit particular stars. 
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1. Introduction 

Mass loss by dust-driven winds of asymptotic giant branch 
(AGB) stars is probably one of the major mechanism 
which recycle material in the Galaxy (e.g. Sedlmayr[l994). 
Most stars (M* < 8M©) will eventually become AGB 
stars and subsequently end their life as white dwarfs sur- 
rounded by planetary nebulae. 

AGB stars are cool (T+ < 3500 K) and luminous (L+ 
of a few 10 3 to a few 10 4 L Q ), and a majority of them 
are pulsating long-period variables (LPVs). The outer lay- 
ers of many AGB stars provide favorable conditions for 
the formation of molecules and dust grains. Dust grains 
play an important role for the heavy mass-loss (up to 
M ~ 10~ 4 M©/yr) of these stars by transferring momen- 
tum from the radiation to the gas. Pulsation causes an 
extended atmosphere where the dust is condensing. The 
dust absorbs the light of the central star and re-radiates 
it at longer wavelengths in the infrared range (A > 2/im). 
The density of material in the circumstellar envelope may 
be so large that the star completely disappears in the vi- 
sual range. 

Send offprint requests to: A.C. Andersen 



Considerable effort has been put into a theoretical de- 
scription of the mass-loss of AGB stars. At the end stages 
of stellar evolution the mass-loss rates become so high 
that they (and not the nuclear burning rates) determine 
the stellar evolutionary faith. The goal is to develop a 
mass-loss description that can be used as input to mod- 
els of stellar evolution and the chemical evolution of the 
Galaxy. 

This paper is the fifth in our series on dust forma- 
tion in winds of long-period variables. The previous mod- 
els presented in Paper I-IV (Dorfi & Hofner 1991 Hofner 
& Dorfi IP9U21 Hofner et al. 119351 Hofner & Dorfi ITM7|> 
are all based on gray radiative transfer while the models 
presented here are calculated using frequency dependent 
radiative transfer for the gas and dust (see Hofner et al. 
US! for details). In this paper we focus on the influence 
of the micro-physical (i.e. optical and chemical) properties 
of dust grains on the winds of AGB stars, to determine to 
what extent the choice of micro-physical parameters ef- 
fects the general mass-loss predictions. 

In Sect. 2 we discuss the hydrodynamical models with 
emphasis on the treatment of the grain formation. Section 
3 describes the different types of amorphous carbon dust 
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used to test the model dependence on the choice of the 
micro-physical parameters. In Sect. 4 we present the re- 
sults and show that the choice of opacity data, the val- 
ues used for the sticking coefficients and the intrinsic dust 
density of the material, all affect the mass-loss rates and 
other wind properties resulting from the models. This has 
an influence on the calculated synthetic colors which are 
compared with observations. Conclusions are presented in 
Sect. 5. 

2. The hydrodynamical model 

Our spherically symmetric hydrodynamical models pre- 
dict the mass-loss rates of AGB stars by treating in de- 
tail the atmosphere and the circumstellar environment 
around pulsating long-period variable stars. This is done 
by solving the coupled system of frequency-dependent ra- 
diation hydrodynamics and time-dependent dust forma- 
tion (cf. Hofner et al. I2002f) employing an implicit nu- 
merical method and an adaptive grid (for details on the 
numerical technique see Dorfi & Feuchtinger 1995). 

2.1. Hydrodynamics and radiative transfer 

In our radiation-hydrodynamical models the stellar and 
circumstellar envelope are described in terms of separate 
conservation laws for the dust, the gas and the radia- 
tion field. The resulting set of nonlinear partial differential 
equations consist of: 

— Equation of continuity (mass conservation). 

— Equation of motion (matter momentum conservation) . 

— Equation of gas internal energy (matter energy conser- 
vation) . 

— 0. moment of the radiative transfer equation (radiation 
energy conservation). 

— 1. moment of the radiative transfer equation (radiation 
momentum conservation) . 

— Moment equations for the dust (grain formation and 
grain growth). 

— Poisson equation (self gravitation). 

To this system of nonlinear partial differential equations 
we add the so-called grid equation which determines the 
locations of the grid points, depending on some critical 
physical quantities to be resolved during the computa- 
tions. 

For the dust component only certain moments of the 
grain size distribution have to be known for a complete 
description of the circumstellar envelope (more details in 
Sect. EHTTJl . 

In contrast to earlier models in this series the models 
presented here describe the radiation field by a frequency- 
dependent treatment of the gas and dust (Hofner 1999 
Hofner et al. I2l)l)2f) based on opacity sampling data of 
molecular opacities (SCAN data base, J0rgensen lf997[) at 
51 frequency points between 0.25 and 12.5 /im. Solving the 
frequency-dependent transfer increases the computation 
time per time step considerably. To keep the computation 



time at a reasonable level, the spatial grid points available 
in the model are reduced from 500 to 100 compared to the 
gray models presented in e.g. Hofner & Dorfi (1997). 

Basically four stellar parameters are needed as input; 
the stellar mass (M*), the effective temperature (T e g), the 
luminosity (£*) of the star and the gas abundance ratio 
{ec/so) °f carbon to oxygen. However, to describe the 
pulsation of the star, two additional parameters are re- 
quired. Pulsation is simulated by a sinusoidal motion of 
the inner boundary Ri n which is located below the stellar 
photosphere. This variable inner boundary is parameter- 
ized by a velocity amplitude (Am) and a pulsation period 
(P) . The luminosity at the inner boundary is variable and 
no mass flux is allowed across this boundary (cf. Hofner 
& Dorfi 1997). As in previous models a perfect gas law 
with 7 — 5/3 and fi — 1.26 is used as equation of state for 
reasons of comparability. 

The dynamical calculations start from a dust free hy- 
drostatic initial model. The radiative pressure on newly 
formed dust initiates an outward motion and the expan- 
sion is followed by the grid out to around 20-30 i?*. At 
this radius the outer boundary is fixed, allowing for out- 
flow. A model evolves for typically 100 years. The model 
calculation is stopped before a significant depletion of the 
mass inside the computational domain occurs. 

2.2. Grain formation in the models 

Dust influences the dynamics and thermodynamics of the 
stellar atmosphere by its opacity. To calculate this opacity 
we need to know the amount of dust present. Dust forma- 
tion proceeds far from equilibrium and it is necessary to 
use a detailed time-dependent description to determine 
the degree of condensation and other relevant properties 
of the grains. 

Dust is formed by a series of chemical reactions in 
which atoms or molecules from the gas phase combine 
to clusters of increasing size. The molecular composition 
of the gas phase determines which atoms and molecules 
are available for the cluster formation and grain growth. 
Dust formation begins with nucleation of critical clusters 
followed by growth to macroscopic dust grains. 

2.2.1. Grain growth 

In our models grain formation is treated by the so-called 
moment method (Gail & Sedlmayr 1988. Gauger et al. 
1990|. The moment method describes the time evolution 
of an ensemble of dust grains of various sizes (includ- 
ing effects due to chemical and thermodynamical non- 
equilibrium) and requires the nucleation rate as external 
input. 

In the moment method the size of a single dust grain 
is expressed in terms of the number of monomers (N) con- 
tained in the grain. A monomer represents the basic ele- 
ment the grain is built of, i.e. a certain species of atoms or 
molecules. The ensemble of the grains is described by the 
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size distribution function f(N,t) representing the number 
densities of dust particles in dependence of their size N. 
The number density f(N,t) of dust grains containing N 
monomers is changed by four processes: creation of grains 
of size N by growth of smaller dust particles and by de- 
struction of larger ones as well as destruction of grains 
with N monomers by growth or evaporation (see Gail & 
Sedlmayr 198||and Gauger et al. 119901 for details). 

As long as the size of the dust grains is small compared 
to the wavelength of the photons in the relevant frequency 
region, which is the case for late-type stars, the optical 
properties of the dust do not depend explicitly on the size 
spectrum f(N,t) but only on a few moments (Kj) of the 
grain size distribution function defined by 



N^ d f(N,t), 



N=N, 



where d denotes the spatial dimension of the grains (d = 3 
for spherical particles and 2 for planar structures) and 2Vj 
is the lower limit of the grain sizes which may be regarded 
as macroscopic in the thermodynamical sense. In the work 
presented here we use the value d = 3 (i.e. assume spher- 
ical grains) and Ni = 1000. 

From the moments Kj it is possible to calculate quanti- 
ties like the total number density of dust grains, the mean 
grain radius and the fraction of condensible material ac- 
tually condensed into grains. 

Considering grains large enough that their thermody- 
namical properties do not depend on the grain size and 
assuming that only molecules with up to a few monomers 
contribute significantly to the growth process the follow- 
ing set of moment equations can be derived (Gauger et al. 
1990) 

dK , (2) 



dt 
dK n 



J 



dt - dr 



+ Np d J 



(!<3<d). 



(3) 



Here Eq. determines the grain production rate. The 
quantity 1 /t is the net growth rate of the dust grains and 
J is the net transition rate per volume from cluster sizes 
N < Ni to N > N[, this can be interpreted as the local 
current density of clusters flowing up-wards in cluster size 
space from the region N < Ni (Gail & Sedlmayr 1988). 

The net growth rate 1/t contains the number densi- 
ties of the chemical species which take part in the dust 
formation process. In our models, the relevant quantities 
are obtained by assuming chemical equilibrium in the gas 
phase at the gas temperature T g 
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where K, denotes the dissociation constant of the molecule 
of the growth reaction and K7 is the dissociation con- 
stant of the molecule involved in the reverse reaction (see 
Gauger et al. I1990|) . f(i) and n ijTO are the number den- 
sities of the i-mers and the molecules containing i-mers 
which contribute to the grain growth, respectively. A\ de- 
notes the (hypothetical) monomer surface area A\ = Aira^ 
where a\ is the monomer radius (see Eq. JSJ), Vth the ther- 
mal velocities of the corresponding growth species and a 
the sticking coefficients. The supersaturation ratio S is 
defined by 



(1) S 



(5) 



which is the ratio of the actual partial pressure of the 
monomers in the gas phase to the vapor saturation pres- 
sure with respect to the dust temperature If the 
thermodynamical conditions allow the formation of dust 
grains the net transition rate J is assumed to be equal 
to the nucleation rate, i.e. the rate at which supercriti- 
cal (stable) clusters are formed out of the gas phase (see 
Sect . I2T21 for details). 

The number densities of the molecules relevant to the 
dust formation are calculated assuming chemical equilib- 
rium between H, H2, C, C2, C2H and C2H2 after the 
fraction of carbon bound in CO has been subtracted. 
Nucleation, growth and destruction of dust grains are sup- 
posed to proceed by reactions involving C, C2, C2H and 
C2H2. The values of the elemental abundances are taken 
to be solar ( Allen ITl)73[) except for the carbon abundance 
which is considered as a free parameter. The dissociation 
constants K,(T) have been extracted from the JANAF ta- 
bles (Stull & Prophet UMJ). 

2.2.2. Grain nucleation 

Nucleation is the first stage of the condensation process 
whereby a vapor transforms to a solid or liquid. This 
phase change requires some degree of supersaturation in 
order to drive the system through the relatively unsta- 
ble reactive intermediates (clusters) between the atomic 
or molecular vapor and the macroscopic solid or liquid 
states. Presently no nucleation rates based on calculations 
of chemical pathways are available for the astrophysical 
problem under consideration here. Consequently, the nu- 
cleation rate which is a function of temperature, density 
and supersaturation (<S) for a particular vapor is often 
calculated by either the classical homogeneous nucleation 
theory (Becker & Doring 1935 Fcder et al. I1966|l or by 
the related scaled homogeneous nucleation theory (Hale 

mm- 

The classical homogeneous nucleation theory was de- 
veloped to describe the nucleation of volatile 1 materials 
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such as water, hydrocarbons or alcohols at relatively low 
levels of supersaturation (S ~ 1.1 — 5.0) and tempera- 
tures (~ 300 K). The theory was as such not developed to 
deal with supersaturated refractory 2 vapors at high tem- 
peratures. The theory describes the formation of critical 
nuclei in a supersaturated vapor by means of thermody- 
namic quantities. The essential basic assumption of this 
approach is that the properties of the clusters in the nu- 
cleation regime are given by the extrapolation of the bulk 
properties even into the domain of very small clusters or 
the interpolation of thermodynamic properties between 
those of the molecules and the solid particles. With these 
assumptions both the thermodynamic functions such as 
entropy and enthalpy and the rate coefficients describing 
cluster formation and destruction become simple analyti- 
cal functions of the cluster size N, which allow a straight- 
forward calculation of the rate of formation of critical clus- 
ters. 

A fundamental result of classical nucleation theory is 
the existence of a bottleneck for particle formation. The 
small unstable clusters which form at random from the 
gas phase have to grow beyond a certain critical size JV* 
which corresponds to a maximum in the Gibbs free energy 
of formation and separates the domain of small unstable 
clusters from the large thermodynamically stable grains. 
The rate of grain formation is determined by the transi- 
tion rate J* between both regions. The existence of such 
a critical cluster size also holds in more realistic theories 
of cluster formation. However, a review of the available 
experimental literature by Nuth & Ferguson (1993) shows 
that no experimental data exists to support the applica- 
tion of classical nucleation theory to the condensation of 
refractory vapors. Refractory vapors seem to condense out 
at different supersaturation ratios than volatile materials. 

The application of classical homogeneous nucleation 
theory in astrophysical discussions of grain formation is 
sharply criticized by Donn & Nuth ( 1985 ) . There was some 
hope that scaled nucleation theory (Hale 11986(1 might be 
a better way to describe the condensation of refractory 
vapors The scaled nucleation theory is a generalization 
of classical nucleation theory by scaling the relevant pa- 
rameters to those of the vapor at the critical temperature 
and pressure and the agreement with experimental data 
for various molecular fluids was rather good (Martinez et 
al. 2001). As a result of this success scaled nucleation the- 
ory was subsequently applied to selected refractory nucle- 
ation data (Hale et al. H989fl with some success. However, 
comparison materials suggest that refractory materials in 
general can not be described as accurately (i.e. significant 
deviations occur for lithium, magnesium and bismuth) as 
the molecular fluids to which the scaled nucleation theory 
was originally applied (Nuth & Ferguson 1993). It was 
shown by Martinez et al. (200TJ that the reason for the 
poor agreement of scaled nucleation theory to certain re- 
fractory materials appear to be, at least in part, the result 
of using an overestimated value for the excess surface en- 

2 A material that vaporises only at high temperature. 



tropy for liquid metals. Martinez et al. l|2001[l conclude 
that refractory materials, as a class, seem to behave dif- 
ferently than the simple fluids studied in the original work 
by Hale ( 1986$ and that the use of bulk liquid properties 
to describe a process involving small metallic clusters is 
problematic and that there therefore is a serious need for 
more and better nucleation data for refractory materials. 

A better description of the nucleation of small refrac- 
tory clusters, which are needed as input for the moment 
method, will most likely have to be guided by experiments. 
For now we will use the classical nucleation theory for 
our model calculations. We are aware that in this way we 
introduce uncertainties to our dust description but with 
no significant improved theory we find it justified to use 
classical nucleation theory as a first crude approximation 
for calculating the nucleation rate Once an improved 
description is derived it will be relatively easy for us to 
change the description of the nucleation rate in the code. 

In the present models the supersaturation ratio S is 
defined as the ratio of the partial pressure of carbon atoms 
in the gas phase divided by the saturated vapor pressure 
of solid carbon (Eq. ©). As shown in Gail & Sedlmayr 
( 1988J S depends on the actual lattice temperature of the 
A^-cluster. 

The value <S = 3 is adopted following Gail & Sedlmayr 
l|1987b|) as the minimum value for grain nucleation to oc- 
cur. As seen from Fig.^ where the supersaturation ratio 
S and the corresponding nucleation rate J* is shown for 
one of the calculated models, the value of S is on the or- 
der of 100 (or at least much larger than unity) in the zone 
where the nucleation rate peaks. 

3. Amorphous carbon dust 

Amorphous carbon particles are considered to be the most 
common type of dust present in circumstellar envelopes of 
carbon-rich AGB stars. The infrared spectra of late-type 
stars generally show a dust emissivity law Q(X) ~ 
with a spectral index of (3 ~ 1 (e.g. Campbell et al.[1976 
Sopka et al. 119551 Martin & Rogers ITW1 Giirtler et al. 
1995). A A" 1 behavior can be expected in a very dis- 
ordered material like amorphous carbon (e.g. Huffmann 
1988} Jager et al. H998|) . Graphite formation in AGB stars 
seem unlikely, because of the absence of the narrow band 
at 11.52 /an in the observed spectra and the overall shape 
of infrared graphite spectra which are proportional to A~ 2 
(e.g. Draine & Lee 11984)) . This is consistent with physi- 
cal considerations, predicting the formation of inhomoge- 
neous grains with crystalline cores surrounded by amor- 
phous mantles (Gail & Sedlmayr I1984|) . 

3.1. Optical properties of dust 

The formation of the dust grains influences the stellar at- 
mosphere in two ways: In the gas phase chemistry, dust 
formation results in a depletion of certain elements, which 
influences the molecular composition of the gas, and con- 
sequently the corresponding opacities. On the other hand, 
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Fig. 1. A snapshot of the supersaturation S and the nu- 
cleation rate J* for model 113djl0pl99. The nucleation 
rate has a sharp maximum around 3i?0 and at this point 
the supersaturation ratio is around 100. 



This means that the opacity of the dust particles requires 
the knowledge of the moment K% and the monomer radius 

trip \ 

~A ' ^ 

47rp grain J 

here A mon is the atomic weight of the monomer (for carbon 
Anon = 12.01115), m p the proton mass and p gra in is the 
intrinsic density of the condensed grain material. 

In an astrophysical context, knowledge of the moments 
Ki for i = 1,2, 3 is needed if the extinction coefficient is 
required and additionally the moments with i = 4, 5, 6 
if both absorption and scattering are required separately 
because the scattering coefficient K sca (A) depends on Kq 
(Gail & Sedlmavr irggljl . 

The photospheric spectral energy distribution of AGB 
stars has its maximum around wavelengths of 1 /im (with 
a sharp decline of the stellar flux toward shorter wave- 
lengths), and both observations and theoretical arguments 
indicate that typical grain sizes in these stars are much 
smaller than 1 /jm. Thus, the limit of particles being small 
compared to the wavelength is valid for a major fraction of 
the spectrum. We therefore do not need to specify details 
about the size distribution of the grains to calculate the 
grain opacities required for model atmospheres. However, 
it may still be necessary to know properties like the size 
distribution and shapes of the grains to compute detailed 
synthetic spectra at very short wavelengths. The complex 
refractive index of the material as a function of wavelength 
can be determined from laboratory measurements. 



dust grains have a rather high mass absorption coefficient 
which often may be comparable to the gas opacity or even 
exceed it. The total opacity of an ensemble of spherical 
dust grains can be formulated as 

/>OC 

k(X) = / a 2 7r Qcxt (a, A) n (a) da, (6) 
Jo 

where n (a) da is the number density of grains in the grain 
radius interval between a and a+da and Q cx t is the extinc- 
tion efficiency, i.e. the ratio for the extinction cross sec- 
tion to the geometrical cross section of the grain (Bohren 
& Huffman 11983(1 . This means that the size distribution 
function of the dust grains has to be known. However, if 
the particle size is small compared to the wavelength of 
the radiation, the small particle limit (SPL), the extinc- 
tion efficiency Q e xt is given by 

QUX) = ^ = 8 i*(^), (7) 
a A \2 + J 

where 9 denotes the imaginary part and m = n+ik is the 
complex refractive index. The dependence of the opacity 
on wavelength and grain size can therefore be separated 
into two independent factors 

/>oo 

«c X P t L = Qcxt(A)^ / a 3 n(a)da = n a 3 Q'^ t (X) K 3 . (8) 
Jo 



3.2. Laboratory measurements of amorphous carbon 

Due to the nature of amorphous carbon dust materials 
they span a broad range of micro-physical properties. 
In the models presented here we have used three differ- 
ent laboratory measurements (referred to as Jdger 1000, 
Jdger 400 and Rouleau, see Table to describe the opac- 
ity of the dust grains formed in the circumstellar envelope. 

The data by Jager et al. I|1998|) are produced by py- 
rolizing cellulose materials at different temperatures. The 
materials are characterize in exemplary detail. In this 
study we have used the data synthesized at 400° C and 
1000° C. The two materials differ in the bonding, where 
the Jdger 400 material has more single bonds (the carbon 
atoms are mainly sp 3 hybridized) while the Jdger 1000 
material has more double bonds (the carbon atoms are 
mainly sp 2 hybridized) . Reflectance spectra were obtained 
of the samples and from these the complex refractive in- 
dex (m) was derived by the Lorentz oscillator method (see 
e.g. Bohren & Huffman HHH Chap. 9). 

Rouleau & Martin (|1991|1 produced synthetic optical 
constants (n and k; m = n + ik) based on measurement of 
sub-micron amorphous carbon particles by Bussoletti et 
al. l(T^S7jl . The particles were produced by striking an arc 
between two amorphous carbon electrodes in a controlled 
Ar atmosphere. 
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Table 1. List of the different dust opacities used. 



Reference 


Material 


Pgr ain 


sp z 


Designation 


Comments 




name 


(g/cm 3 ) 


% 


in this paper 




Jager et al. ©98 ) 


cel400 


1.435 


67 


Jager 400 


completely amorphous 


Jager et al. 


cellOOO 


1.988 


80 


Jdger 1000 


contains graphite (2 nm) crystallites 


Rouleau & Martin (1991) 


AC2 


1.85 




Rouleau 





3.3. Grain equilibrium temperature 

The presence of dust grains influences both the momen- 
tum and the energy balance of the atmosphere. We as- 
sume complete momentum coupling of gas and dust, which 
means that the momentum gained by the dust from the ra- 
diation held is directly transferred to the gas. On the other 
hand, the transfer of internal energy between gas and dust 
is negligible compared to the interaction of each compo- 
nent with the radiation field (cf. Gauger et al. I1990[) . We 
therefore assume that the grain temperature is given by 
the condition of radiative equilibrium 



X.i J - Xs S d = 



Ta 



xj_\ 

XsJ 



(10) 



where the source function is equal to the Planck function 
Sd = B(T d ). The frequency- integrated dust opacities are 
defined by 



X.J = 



Xs 



J v Xv J_v dv 
J J„ dv 

S v Xv S v dv _ j jL Xu B u (T d )dv 

J Sydv f B u (T d ) dv 

j v Xv B v {T d )dv 
B(T d ) 



(11) 
(12) 



where J v is the radiation energy density, and the radiation 
temperature T r = (J^/cr) 1 / 4 . 

For a gray opacity we would obtain \J = Xs for the 
opacities and the temperatures T d = T r as in previous 
models. Using frequency-dependent radiative transfer the 
opacities x.J an d Xs differ. For the data shown in Fig.|2 
the absorption coefficients decrease with increasing wave- 
length in the infrared region of the spectrum, leading to 
Td > T r for all data sets shown. The difference between Td 
and T r becomes larger with an increasing slope of Q' QXt W 
and the difference may reach several 10 2 K in the dust 
formation zone (Hofner et al. 2002). In other words, the 
steeper the dependence on wavelength, the larger the dif- 
ference between the equilibrium grain temperature and 
the radiation temperature. Therefore, the steeper slope of 
Jdger 400 data results in a relatively high grain tempera- 
ture compared to Jdger 1000 and Rouleau. 

4. Results 

One of the main reasons for the huge mass-loss of AGB 
stars seems to be the presence of newly formed dust grains. 
The strong shock waves in the stellar atmosphere cause 




Jager 10001 
Jager 400 : 
Rouleau - s 



io°L 

0.1 



1.0 10.0 100.0 

Wavelength (nm) 



1000.0 



Fig. 2. Wavelength dependence of the dust absorption ef- 
ficiency Q' C xt (see Tabled for dust annotation). 



a levitation of the outer layers. The cool and relatively 
dense environment which results from the levitation pro- 
vides favorable conditions for the formation of molecules 
and grains. Due to its high opacity and the resulting ra- 
diative pressure, the dust has a strong influence on the 
structure of the atmosphere and the wind properties. 

4.1. Influence of the dust extinction coefficient 

The influence of the dust extinction coefficient on the 
winds of the dynamical models is shown in Tablc|5] and 
Fig-El The models depend on the choice of laboratory 
measurements of Q' cxt {\) as already shown for gray mod- 
els in Andersen et al. ( 1999). In contrast to the gray mod- 
els, for the frequency-dependent models both the absolute 
value and the slope of the dust opacity data as a func- 
tion of A become relevant. The absolute value of the grain 
opacities mainly affects the terminal velocity of the winds, 
while the slope has significant influence on the grain tem- 
perature as discussed in Sect . 13.31 In the models where the 
Jdger 400 data is used for the opacity of the grains, the 
slope of the extinction efficiency as a function of A dictates 
a high grain temperature which prevents dust formation 
for the chosen stellar parameters. Models with identical 
stellar parameters using Jdger 1000 and Rouleau data on 
the other hand develop dust-driven winds. 

Figure [3] and Table show that it is not possible to 
distinguish from near infrared photometry alone between 
the different models, since the hotter model using the 
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Fig. 3. The colors (J— H) vs. (H— K) of the four models 
113djl0pl99 (o), 110djl0pl99 (A), HSdrouplSS (x) and 
110droupl85 (*). There are two differences between the 
two sets of models. (1) The assumed value of the luminos- 
ity (Z/*) and corresponding effective temperature (T*), and 
(2) the dust opacity data used (Jager 1000 and Rouleau). 
All other model values are kept constant. See Table for 
details. 

Jager 1000 opacity data has similar colors as the cooler 
model using the optical properties from the Rouleau data. 
This is a consequence of the fact that these two mod- 
els have circumstellar envelopes with comparable optical 
depth at these wavelengths as can be demonstrated by 
the following rough estimate; the optical depth of a given 
layer dr can be expressed in terms of the extinction coef- 
ficcnt, the local degree of condensation / c , mass- loss rate 
and flow velocity by 

dr = ndr = ir afQ' cxt (X) K 3 (r) dr (13) 
« a i Q'cxt fc Pgas dr 

K a i Q'cxt fc — -4 dr 
u r z 

where we have used Eq. JSJ, K% oc f c p gas and p gas = 
M /Anr 2 u to replace the gas density p gas - To estimate the 
relative optical depths of the circumstellar envelopes of 
different models we assume that we can replace f c M/u 
by the the time-averaged values (f c )M/(u) at the outer 
boundary. For models 113djl0pl99 and 110droupl85 the 
differences in mass-loss rate, degree of condensation and 
outflow velocity as well as the extinction efficiency and 
the monomer radius compensate each other in such a way 
that the quantity afQ' ext {f c )M / (u) is comparable for both 
models while the value for 110djf0pl99 differs by about 
a factor of two. Although it is impossible to distinguish 
the models If 3djI0pI99 and 110droupI85 by their near in- 
frared colors alone it would be possible to do this by a com- 
bination of photometry and high-resolution spectroscopy 
due to their significantly different outflow velocities (of 
about a factor of four). 



Fig. 4. The colors (J— H) vs. (H— K) for the four mod- 
els 1^10^199 (o), 113djl0p225 (*), lI3droupl85 (x) and 
H3droup225 (A). We see the results for two different opac- 
ity data (Jager 1000 and Rouleau). Two different models 
for each opacity data are shown. The difference between 
the two models with the same opacity is the assumed value 
of the intrinsic density (p gra in) of the dust material formed 
in the model. The values correspond to the one deter- 
mined in the laboratory for the material and the value for 
graphite (2.25 g/cm 3 ) which was used in previous models. 
See Table 13 for details. 

4.2. Influence of the intrinsic dust density 

When the moment method was developed by Gail & 
Sedlmayr l|I988[l the intrinsic density for the amorphous 
carbon material used in that generation of models (Maron 
1990J) was not known. The value of graphite (2.25 g/cm 3 ) 
was therefore assumed. This value was later used in 
most existing models based on the moment method (e.g. 
Fleischer et al. lTM?! Hofner & Dorfi lTM?)) . 

In the previous section we have described models which 
use optical properties represented by three different amor- 
phous carbon materials for which the intrinsic densities 
have been measured in the laboratory. Here, we compare 
them with three models where the same optical proper- 
ties were used but where we assumed the higher value for 
the intrinsic density of graphite instead of the measured 
value of the material (while keeping all other parameters 
constant) . 

The result of using the higher density of graphite in- 
stead of the correct values can be seen in Table |3 and 
Fig. 0| Figure0]demonstrate that the models become much 
redder when the intrinsic density of graphite is used in- 
stead of the respective values for amorphous carbon. The 
large increase in the mean degree of condensation (/ c ) 
at the outer boundary implies that more dust is formed. 
This results in an increased radiation pressure on the dust 
grains and the models therefore show a higher mean veloc- 
ity (u) at the outer boundary. At the same time the mass- 
loss rates increases so that a\(f c )M/(u), which is a mea- 
sure of the optical depth of the envelope, is significantly 
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Table 2. Influence of the extinction efficiency of the dust. Model parameters: M* = 1.0 Mq, sticking coefficient 
«c = 0.37, cec 2 — 0.34, ac 2 H = 0.34 and ac 2 n 2 — 0.34, ec/so = 1-4, period P — 650d, piston velocity Am p = 4km/s, 
luminosity L*, temperature T*, radius dust opacity data Kdust, intrinsic dust density p gra in; Results: Mass loss rate 
M, mean velocity at the outer boundary (u), mean degree of condensation at the outer boundary (/ c )- 
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0.10 
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1.99 
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110droupl85 
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2.3 ■ 10 -6 
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Jager 400 


1.44 











Table 3. Influence of the intrinsic density of the dust. Model -parameters: Same as TableEl 
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Table 4. Influence of the sticking coefficient for the dust formation. Model parameters: Same as Table|2 except for 
the sticking coefficients: «c, ac 2 , ac 2 H and ac 2 H 2 , dust opacity Kdust = Rouleau, and the intrinsic dust density 
Pgrain = 1-85 g/cm 3 for all models. 
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[M Q /yr] 
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increased the monomer radius a\ decreases (see Eq. (0). 
This influences both the grain growth and the dust opac- 
ity. The net growth rate 1 /t and the opacity k dependent 
on the monomer radius as 1/r tx a\ (through the monomer 
surface A{) and k cx af , respectively. A reduction of a\ de- 
creases the dust mass absorption coefficient stronger than 
the growth rate, leading both to a slower acceleration and 
a more efficient growth of the grains. The opacity both 
depends on a\ (which decrease slightly) and the degree 
of condensation / c (which increases) with the later effect 
dominating, which leads to a higher final outflow velocity 
(u) of the wind, despite a reduction of the dust absorption 
per gram of dust. 

4.3. Influence of the estimated sticking coefficient 

The sticking coefficient (also called the reaction efficiency 
factor) a enters into the net growth rate of the dust grains 



higher for both models using the density of graphite com- 
pared to the respective models with the consistent dust 
densities. 

Even a small increase of about 10% in the density of 
the dust material (as it is the case from model 113djlO/9l99 
to 113djl0p225) results in a doubling of the degree of con- 
densation and a substantial increase of the outflow velocity 
(u) and the mass-loss rate. For the models 113droupl85 
and 113droup225 where the difference in the value used 
for the intrinsic dust density is about 20%, the estimated 
mass-loss rates differ by almost a factor of two. For the 
models using the Jager 400 material (113dj04pl44 and 
U3dj04p225), where the difference is almost 40%, using 
the measured material value instead of the higher value 
for graphite results in a model that will not develop a 
wind at all. 

The intrinsic dust density of the dust material /9 gra i n is 
so significant for the obtained results because when it is 
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Fig. 5. The colors (J— H) vs. (H— K) for the two models 
liadroupISS (*) and 113droupl85al0 (A). The only dif- 
ference between the two models is the assumed value of 
the sticking coefficient. The larger the sticking coefficient 
is assumed to be, the more efficient will the grain forma- 
tion be and as a consequence the models will look redder 
when observed in (H-L), (J-H) or (K-L). See Table H 
for details. 



(Eq. (JSJ). However, a is not definitely known as long as we 
do not know explicitly the sequence of chemical reactions 
responsible for the dust formation. To demonstrate how 
uncertainties in the value of a will influence the results 
of the models we have varied this parameter in otherwise 
identical models. 

In the case of carbon dust formation the most impor- 
tant growth species is expected to be C2H2, but also C, 
C2 and C2H contribute to the growth (Gail & Sedlmayr 
1988). In our earlier models we have used the values 
ac = 0.37, «c 2 = 0.34, qc 2 h = 0.34 and ac 2 H 2 — 0.34 to 
describe the grain growth. 

Gail & Sedlmayr Q1984[) used the value a c = 0.3 
adopted from Landolt-Bornstein (1968). Later Gail & 
Sedlmayr (1988) argued that the sticking coefficient a 
must be on the order of unity, because it is expected that 
neutral radical reactions play a dominant role in the for- 
mation process of carbon grains. However, if substantial 
energy barriers are involved in the reactions a may well 
be less by several orders of magnitude. It is important 
to remember that if a is small the assumption that nucle- 
ation can be treated as a time-independent process will no 
longer hold. Small values of a should be accompanied by 
a time-dependent treatment of the dust nucleation (Gail 
& Sedlmayr \MB- 

Salpeter H1973[l has shown that the latent heat released 
when a monomer attaches itself to an 7V-mer may lead to 
a small sticking probability if ./V is small. Thus the sticking 
coefficient a for clusters with N ss N* may be considerably 
smaller than the sticking coefficient for clusters with N 3> 
N*. The sticking coefficient is of the order unity for bulk 
material (Found lrbT72j) . 



To test the assumptions for the sticking coefficient 
we have calculated four models for a star with = 
13000 L , T* = 2700 K and with the sticking coefficient 
(a) varying from 0.2 — 1.0 and two cooler models with 
L* = 10000 L Q and T* = 2600 K. 

It is clear from Tableland Fig.[5](showing how the col- 
ors of two models depend on the chosen value for the stick- 
ing coefficient) that using the value of 1 as suggested by 
Gail & Sedlmayr (1988 ) compared to our previous values 
of 0.34 — 0.37 results in a noticeable reddening of the colors 
of the stellar model. Fig.[|J] shows the mean outflow veloc- 
ity, the gas density, the gas temperature and the mean 
degree of condensation for four different phases of the two 
different models 113droupl85a and 113droupl85al0. The 
higher value of the sticking coefficient increases both the 
degree of condensation and the outflow velocity by about 
a factor of two. In addition the mass-loss rate increases, 
leading to a higher optical depth of the circumstellar en- 
velope. 

From Fig.HJlit can be seen that dust formation occurs 
beyond « 2R+ independent of the value assumed for the 
sticking coefficient since the onset of condensation is de- 
termined mainly by the temperature. With the low choice 
of sticking coefficient only about 10% of the condensable 
carbon material present in the gas actually condenses into 
grains, while for the much more favorable choice of stick- 
ing coefficient (a = 1) twice as much material condenses 
into grains. A complete condensation of carbon grains is 
prevented by the rapid velocity increase after the onset 
of avalanche nucleation and grain growth, and the subse- 
quent rapid dilution of the gas. 

4.4. Influence of the surface tension for the dust 

In classical homogeneous nucleation theory, the surface 
tension of the grain material is used to describe the gain 
of enthalpy by forming a grain out of N monomers in the 
gas phase. 

The surface tension Ograin °f amorphous carbon is not 
known from laboratory experiments, it has therefore be- 
come custom to use values for graphite, e.g. one of the 
values given by Tabak et al. I|1975[) . One problem however 
is that there are huge variations in the values determined 
due to the anisotropy of graphite. 

As already demonstrated by Tabak et al. I|1975f) . vary- 
ing the value of the surface tension may produce an enor- 
mous change in the nucleation rate. To determine how 
significant the prescribed value for the surface tension 
is for our results we have varied it around the value of 
1400 erg/cm 2 which was used in previous models (see 
Table0). 

For the particular model, altering the value of the sur- 
face tension of the dust grains by 28% around the value 
of graphite will make the difference between obtaining 
mass- loss or not. The value of the surface tension also 
has a substantial influence on how much of the available 
material in the circumstellar envelope will condense into 
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Fig. 6. Models 113droupl85 (left) and 113droupl85al0 (right). For four different phases of the two models (a) the 
outflow velocity of the gas, (b) the gas density, (c) the gas temperature and (d) the degree of condensation is shown. 
The only difference between the two models is the value used for the sticking coefficient. It is evident that the choice 
of sticking coefficient has a significant influence on the predicted outflow velocity and the degree of condensation. 
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Table 5. Influence of the surface tension for the dust for- 
mation. Model parameters: Same as for model llSdrou/3185 
in Tabled] but with different values for the surface tension 



Model 


C grain 


M 


<«> 


</c) 




[erg/cm 2 ] 


[M /yr] 


[km/s] 




113droupl85al0 


1000 


8.5 ■ 10~ b 


35 


0.76 


113droupl85 


1400 


4.9 ■ 10" 6 


7.4 


0.10 


113droupl85al8 


1800 









1.8 P 




J-H 

Fig. 7. The colors (J - H) vs. (H - K) for the three mod- 
els 113djl0pl99 (A), 113droupl85al0 (□) and 110djl0pl99 
(O)j an d three comparable observed carbon Miras R Vol 
(+), R For (*), R Lep (x) (observations from Whitelock 
et al. lTOi . 

dust grains. For comparison, the measured surface tension 
for other materials are: Fe (er gra in = 1400 erg/cm 2 ), MgS 
(ograin = 800 erg/cm 2 ) and SiO (cr gra in = 500 erg/cm 2 ), 
see Gail & Sedlmayr 

4.5. Comparison with observations 

Near-infrared colors including JHKL 3 magnitudes were 
calculated for a selection of the models. The filter zero 
points were calculated from a Vega model of Dreiling & 
Bell (1980), under the assumption that the Vega model 
has 0.0 mag in all filters. 

Whitelock et al. (1997) present JHKL light curves for 
11 large-amplitude carbon variables. For a comparison be- 
tween models and observations we picked the carbon Miras 
R For, R Lep and R Vol. These stars have moderately 
thick dust shells and colors comparable with our models 
(see Fig.[7Jl . According to Olofsson et al. C5M| these stars 
have outflow velocities of about 16 — 19 km/s. From our 
set of models we selected those with comparable velocities 
(see Table EJ). We preferred outflow velocities over other 

3 Central filter wavelength [/im]: J = 1.22, H = 1.63, 
K = 2.19, L = 3.45. 



wind properties when selecting the models for this com- 
parison, since this is the most directly observable quantity. 
The colors on the other hand contain entangled informa- 
tion about the mass-loss rates, the ouflow velocities, the 
degree of condensation, and the optical properties of the 
grains. 

We have not selected the stars for their derived stellar 
parameters but for their observable wind properties. The 
observed stars (R For, R Lep and R Vol) and the models 
(113droupl85al0 and 110djlO j ol99) have similar dust mass 
loss rates and near infrared colors while the gas mass-loss 
rates for the observed stars and the models differ signifi- 
cantly (see TableEJ). We preferred the dust mass- loss rates 
over the gas mass-loss as a criterion for the comparison 
since the near IR colors are strongly affected by the dust 
opacity, as demonstrated by the figures in the preceeding 
section. 

The dust mass-loss rates of the observed stars were 
estimated using the 60 fim IRAS flux. In general there is 
expected to be a relation between the dust mass-loss rate 
and the strength of the dust emission. Following Sopka et 
al. (1985) Olofsson et al. (1993) derive the dust mass-loss 
rates of these stars using 

M dust = 2.2 x 10- 15 S 60 ^ 2 £~ - 5 M yr-\ (14) 

where Sqq is the IRAS 60^m flux in Janskys, Vd is the 
dust expansion velocity in km/s, D is the distance to the 
source in parsecs and L is the luminosity in units of Lq . 
For the dust expansion velocity Olofsson et al. (1993} use 
two different values, one where the upper limit to the drift 
velocity is obtained by considering the radiation force on 
the grains and the drag force due to gas-grain collisions 
in the limit of supersonic motion. The lower limit for the 
dust expansion velocity is derived by neglecting drift. As 
a consequence two different dust mass-loss rates are given 
by Olofsson et al. I|1993[l . one for the upper limit M^^ r 
and one for the lower limit M l d °^ t er . 

The gas mass-loss rates of the observed stars given 
in Table HO were derived by Olofsson et al. ( 1993) based 
on observed CO emission using equation (5) of Knapp & 
Morris l|1985|) which has the form 

where n-He/H is the total number density of 
Helium/Hydrogen, T m b is the main-beam brightness 
temperature, v e represent the gas expansion velocities, D 
is the distance to the source, A(B, J) is a quantity that 
depends on the beam size used, B, and the transition 
observed, J — > J — 1, and fco is the abundance of 
CO with respect to H 2 . Schoier & Olofsson l|2001|) find 
that the gas mass-loss rates resulting from Eq. (|15|l are 
systematicly underestimated compared to a more detailed 
radiative transfer analysis (see their Fig. 8), however, for 
these particular stars the values agree resonable well. 

In the models discussed here we obtain higher gas-to- 
dust ratios than the values given by Olofsson et al. 0993). 
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Table 6. Data for three observed carbon-rich Mira's with similar colors and outflow velocities as the calculated models 
(Tables and EJ. The observed velocities and mass- loss rates are from Olofsson et al. {1993 ). 
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This could be a consequence of the fact that some of the 
stellar parameters differ from these of the models. As men- 
tioned above, the models and stars for the comparison 
were chosen for similar wind velocities and dust mass-loss 
rates, not for comparable stellar parameters. In this sense, 
we are rather comparing the observable properties of the 
circumstellar dust shells than of the stars and their sur- 
rounding gas envelopes. Therefore, the comparison shown 
here should not be considered as an attempt to fit obser- 
vations of individual stars but to investigate certain wind 
properties of our models. 

While the average colors of the models are compara- 
ble to the average colors of the observed stars, see Fig.0 
the temporal variation of the colors for a given model 
is smaller than the variations observed in the individual 
stars. This is probably a consequence of relatively small 
bolometric luminosity variations in the models. 

Figure shows a comparison of models and observed 
stars for the colors J-H vs. H-K. A qualitatively similar 
relation between models and observations is obtained for 
the colors K-L vs. H-K. 

5. Conclusions 

We have investigated in detail how the predicted wind 
properties of carbon-rich AGB stars are influenced by the 
choice of micro-physical dust parameters i.e. the optical 
properties of the dust, the intrinsic dust density, the as- 
sumed sticking coefficients and the surface tension of the 
grain material (these two last parameters control the effi- 
ciency of the dust formation). 

For the theoretical predictions of mass-loss it is impor- 
tant to know how the uncertainty in the chosen dust pa- 
rameters affects the obtained results. Varying the micro- 
physical parameters within the range typical of possible 
materials can change the value for the mean outflow ve- 
locity of the gas and dust as well as the predicted degree 
of dust condensation by a factor of ten and the predicted 
mass-loss by a factor of four. In the transition region be- 
tween models with and without mass-loss the choice of 
micro-physical parameters is vital for whether a particu- 



lar set of stellar parameters will give rise to a dust-driven 
mass-loss or not. 

The main source of momentum for the stellar wind is 
the radiation pressure on dust. The radiation pressure on 
the dust and the radiative equilibrium grain temperature 
is determined by the wavelength dependence of the grain 
extinction efficiency. The steeper the dependence of wave- 
length, the larger the difference between the equilibrium 
grain temperature and the radiation temperature. The ra- 
diation pressure on the other hand is proportional to the 
flux mean opacity which both depends on the slope of the 
extinction efficiency as a function of wavelength and on its 
absolute value. The latter may differ by almost an order 
of magnitude for different types of amorphous carbon in 
the critical region around 1 /zm. The density of the grain 
material has to be chosen consistently with the grain ex- 
tinction efficiency. 

The surface tension of the grain material and the stick- 
ing coefficients are very significant for the calculated rates 
at which grains are formed out of the gas (nucleation) and 
at which new material is added to existing grains (grain 
growth). Even a moderate variation of the values within 
the range expected for possible materials has noticeable 
consequences for the properties of the dust-driven stellar 
winds, including the resulting near-infrared colors. 

The colors are similar to stars with comparable dust 
mass-loss rates and outflow velocities, however, we have 
not attempted to fit any individual stars. 
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